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The electrodynamical response of the edge of a compressible Quantum Hall system affects tunnel- 
ing into the edge. Using the composite Fermi liquid theory, we derive an effective action for the edge 
modes interacting with tunneling charge. This action generalizes the chiral Luttinger liquid theory 
of the Quantum Hall edge to compressible systems in which transport is characterized by a finite 
Hall angle. In addition to the standard terms, the action contains a dissipative term. The tunneling 
exponent is calculated as a function of the filling fraction for several models, including screened 
and unscreened long-range Coulomb interaction, as well as a short-range interaction. We find that 
tunneling exponents are robust and to a large extent insensitive to the particular model. We discuss 
recent tunneling measurements in the overgrown cleaved edge systems, and demonstrate that the 
profile of charge density near the edge is very sensitive to the parameters of the system. In general, 
the density is nonmonotonic, and can deviate from the bulk value by up to 30%. Implications for 
the correspondence to the chiral Luttinger edge theories are discussed. 
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I. INTRODUCTION 



The edge of a Quantum Hall (QH) system attracts a lot 
of interest because it provides an example of a one dimen- 
sional non-Fermi-liquid. The theoretical picture of QH 
edge was first developed for odd-denominator Landau- 
level filling fractions v that correspond to incompressible 
QH states Q| . It involves one or several interacting chiral 
Luttinger liquid modes. The most prominent feature of 
Luttinger liquid is the power law character of the Green's 
function. 

A power like Green's function leads to a power law 
in the tunneling-current -voltage dependence: / ~ V a . 
The tunneling exponent a has been extensively studied 
theoretically for the principal filling fractions of Laugh- 
lin and Jain hierarchies pl[|. For Laughlin states with 
v =1/ (2k + 1) the edge is described by one chiral mode 
and tunneling current I ~ V 2k+1 is predicted fO]. The- 
ories of the edge with v ^ 1/ (2k + 1) involve more than 
one mode. In the multi-mode case the results are quali- 
tatively different for the modes going all in one direction, 
and for modes going in the opposite directions. 

For comoving edge modes, the tunneling exponent is 
universal and does not depend on the character of inter- 
action between the modes. For example, this is the case 
at the Jain filling fractions v — n/(np+ 1) with positive 
integer n and even p, where Wen jjj finds / ~ V p+1 . On 
the other hand, for the edge described by modes going in 
opposite directions, the tunneling exponent depends on 
the interaction strength. In this case, it is also important 
to take into account the effects of disorder Q . The point 
is that relaxation between the modes due to scattering by 
disorder mixes the modes, and at sufficiently high disor- 
der effectively forms a single charged mode with universal 
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tunneling exponent. For example, for the Jain fractions 
n/(np+ 1) with n < and even p > 0, Kane, Fisher, and 
Polchinski § found I ~ yp+i-2/M. 

Tunneling experiments probing the physics of the QH 
edges were first attempted using conventional split gate 
devices ||], after which a new generation of 2D systems 
was developed by using the cleaved edge overgrowth tech- 
nique |Q. In these structures it is possible to study tun- 
neling into the edge of a 2D electron gas from a 3D doped 
region. In this system one can create a 2DEG with a very 
sharp edge, with residual roughness of an atomic scale. 
High quality of the cleaved edge system permits to ex- 
plore tunneling in both incompressible and compressible 
QH states §§. 

In the first experiment Q, for v = 1/3 it was found 
that the tunneling conductivity is non-ohmic, / ~ V a , 
with the exponent a ~ 2.7, quite close to the theoretical 
prediction a = 3. After that, it was observed that 
the power law / ~ V a holds for both incompressible and 
compressible QH densities, in the range 0.25 < v < 1. 
The exponent a was found to be reasonably accurately 
given by a simple formula: a = 1/v. Interestingly, this 
dependence did not agree with the predictions of chiral 
Luttingcr models, except for a special point v = 1/3. 
Moreover, it was quite surprising that the power law 
is equally well obeyed by both compressible and incom- 
pressible values of v. 

These findings prompted interest in the problem of 
tunneling into the edge of a compressible QH system. 
A good description of the compressible QH states is pro- 
vided by the composite fermion theory. This theory [Q is 
constructed at v = 1/2 and other rational v with even de- 
nominator, and is used to map the problem of fractional 
QH effect onto the integer QH problem for new quasi- 
particles, composite fermions interacting via statistical 
Chern-Simons gauge field. In the composite fermion pic- 
ture, an electron is described as a fermion carrying vor- 
ticity represented by a quantized gauge field vortex. For 
densities close to the half-filled Landau level the vortex 
has p = 2 flux quanta. The theory of composite fermions 
withp = 2 describes the interval of densities 1/3 < v < 1. 
At smaller densities 1/5 < v < 1/3 composite fermions 
with p = 4 are used, etc. 

The theory of tunneling into a compressible QH edge 
|Tof which uses composite fermions to describe the QH 
system predicts the power law / ~ V a , with a being a 
continuous function of v. 
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where p xy — (h/e 2 )v 1 is the Hall resistance of the 

2DEG and p x °y = (/i/e 2 )(z/ _1 — p) is the Hall resistance of 
composite fermions moving in an effective magnetic field 
B c ff = B — p^qu, where n is electron density. The result 
(||) describes the system in the limit p xx — ► 0. The de- 
pendence of a on 1/f is monotonic, and is characterized 
by plateaus in the intervals 2 < 1/v < 3, 4 < 1/^ < 5, 



etc. (see Fig. || below). The plateaus are connected by 
straight lines with slope 2. The cusp- like singularities 
predicted in the dependence a(v) at v = 1/2, 1/4, etc., 
are somewhat smeared when p xx is finite (see Eq. (|6§1 ) 
and Fig. |). 

Interestingly, these results match exactly the predic- 
tions of the chiral Luttinger liquid theory for Jain series 
of incompressible states. Although formally this theory 
lacks continuity in the filling fraction, starting from a 
new set of edge modes for each given filling fraction, the 
exponents a(v) obtained by Wen for v = n/(pn + 1) and 
by by Kane, Fisher and Polchinski for v — n/(pn — 1) 
fall on the continuous curve (Q) found in the composite 
fermion calculation. The exponents of Wen fall on the 
plateaus, while those of Kane et al. fall on the straight 
lines connecting the plateaus. 

However, the disagreement with the experimentally 
measured a(v) requires an explanation. Recently a num- 
ber of theories were proposed trying to resolve this issue 
In one approach, described by Conti and Vignale [11 
Han and Thouless |l2| and Zulicke and MacDonald [13 
tunneling is studied by using a hydrodynamical theory of 
a compressible QH edge, in which the nature of under- 
lying quasiparticles is essentially ignored. From such a 
treatment the desired relation a = 1/ v em erges readily, 



as we will discuss in detail below in Sec. 
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and at the end 

of Sec. III. However, this approach ignores the contribu- 
tion to the electron Green's function of the quasiparticles 
in the QH state, and thus it is in contradiction with the 
presently existing microscopic picture of the QH effect. 

Another line of thought, developed by Alekseev et al. 
|]l5| , is that the experimental system is not what it is 
assumed to be. In particular, it is proposed that in- 
stead of a clean edge the real system contains many lo- 
calized states in sufficient proximity to the edge. Then, 
if one assumes that the tunneling rate bottleneck cor- 
responds to tunneling from the doped region into a lo- 
calized state, and that the density of localized states is 
sharply peaked about the Fermi energy, one finds the de- 
sired result a = \/v. The reason is that in the problem 
involving a localized state no conversion of an electron 
into a quasiparticle is required, and the only effect to be 
considered is a shakeup of the edge plasmon mode, the 
effect equivalent to the x-ray edge problem in the Fermi 
liquid. However, it is not clear why the density of lo- 
calized states should be peaked at the Fermi energy in 
the actual samples. An apparently similar idea has been 
developed earlier by Pruisken et al. jl4| using quite elab- 
orate methods which we have not been able to follow in 
detail. 

Also, a theory using composite fermions was proposed 
by D.-H. Lee and X.-G. Wen [|~7j in which both charged 
(edge plasmon) and neutral (quasiparticle) modes are in- 
cluded. It was assumed, however, that the velocity of 
the charged mode is much larger than the neutral mode 
velocity. In this case, there exists an intermediate energy 
regime in which only the charge mode dynamics is im- 
portant, while the neutral mode response can be ignored. 
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In this energy domain one obtains a = 1/v. It should be 
pointed out, however, that the ratio of the charged and 
neutral mode velocities is of order of In w/r s , where w is 
the distance from the edge to the doped region, and r s is 
the screening radius. Optimistically, the ratio w/r s can 
be as large as 10 which is still not enough to explain the 
power law demonstrated in a wide range of 2.5 decades 
in the bias voltage. 

Another approach trying to rationalize the mea- 
sured tunneling exponent a ~ 1/v was proposed by 
Khveshchenko jL8|. This theory is based on composite 
fermions and is similar in its assumptions to Rcf. PQ| and 
to the present work. However, the calculated tunneling 
exponent is 1/v up to a frequency dependent logarithmic 
correction small in p xx /p xy . We believe that this is due 
to an inconsistency of the analysis ignoring important ef- 
fects accounting for dynamics of free composite fermions. 
One can see that by comparing Eq.7 of Ref. |Q with our 
Eq.^J, and noting that the term describing the free com- 
posite fermion response is missing in Ref. |18|. 

In addition to this controversy, the theory by Lopez 
and Fradkin [|l6j seems to abandon the entire theoretical 
picture of the multi-mode QH edges proposed in Refs. 

ffor the incompressible Jain fractions. Unlike Ref. 
the authors of Ref. Jl6| are not using a microscopic 
mechanism for eliminating the neutral propagating edge 
modes. The construction proposed in Ref. (lq ] involves 
only one charged mode plus two auxiliary Klein factors 
which do not constitute additional propagating degrees 
of freedom. In that, the approach of Ref. |L6|] can be 
comapred to the conventional quantum Hall edge theo- 
ries jjJU in which the velocity of neutral mode is exactly 
zero. If true, this would lead to the a — 1/v dependence 
at arbitrarily low energies. However, it is presently un- 
clear whether the picture of the neutral mode with zero 
velocity can be justified microscopically. 

What complicates the controversy even further is the 
recently presented evidence of a plateau-like feature ex- 
hibited by a(v) in some cleaved edge samples j[s|. The 
value of v near which the dependence a(v) flattens out is 
however quite close to 1/3, whereas the expected plateau 
interval is 2 < v^ 1 < 3. This discrepancy may be ex- 
plained by solution of the electrostatic problem near the 
edge (see Sec. [v] below and Ref. jl9|]) which shows that in 
a wide region adjacent to the edge the density exceeds the 
bulk value by about 20 — 30%. Because of this behavior 
of the density profile, the feature in a(v) observed near 
^buik = 1/3 may correspond to somewhat higher density 
near the edge, with v somewhere between 1/3 and 1/2. 

One other complication is that the analysis of the elec- 
trostatic problem shows that the density profile near the 
edge can be nonmonotonic and, in general, depends quite 
sensitively on the system parameters. This observation 
can make the relation with the theories assuming con- 
stant filling factor somewhat indirect. At present, the 
matter is obviously far from being resolved, and more 
experimental and theoretical studies are needed to clar- 
ify the situation. With this in mind, in this article we 



present an alternative derivation of the results obtained 
in Ref. ficfl , demonstrating their robustness and estab- 
lishing a more direct connection with the chiral Luttinger 
theories of the QH edge. 

The basis of our analysis will be the theory of com- 
posite fermions [Q] . We assume that noninteracting com- 



posite fermions are characterized by pi° x and p x u y ' which 
may depend, e.g., on the filling fraction. The measured 

resistivities are then p xy = p x ° y + ph/e 2 and p xx = p x ° x , 
where p is the number of flux quanta attached to an elec- 
tron (p = 2 for 1/3 < v < 1). The tunneling current 
is expressed in a standard way in terms of the electron 
Green's function. We derive the relation between Green's 
functions of an electron and of a composite fermion, and 
compute the former using a "factorization approxima- 
tion." In this analysis the effects of shaking up slow 
electromagnetic and Chern-Simons gauge field modes are 
separated out. As a result, the tunneling current is ex- 
pressed in terms of electromagnetic response functions 
and the number of flux quanta p. The theory predicts a 
power law / ~ V a with a continuous dependence of the 
tunneling exponent a on the filling fraction. As far as 
tunneling into the edge is concerned, there is no qualita- 
tive difference between compressible and incompressible 
states. The "Luttinger liquid-like" behavior in the edge 
tunneling emerges when the Hall angle is close to tv/2, for 
both compressible and incompressible electron systems. 

The paper is organized in the following way. In Sec. [ B 
we review the approach of Ref. based on a semiclassi- 
cal phase factor analysis of the Green's function. This is 
done with the purpose of motivating and providing con- 
nection with the subsequent discussion of the effective 
action formalism. In Sec. |n| we begin laying out the ba- 
sic approach of the present theory of tunneling. At low 
energy, the most important effect is the shake-up of long 
wavelength modes corresponding to spreading of the tun- 
neling charge. To describe it, one can use a semiclassical 
method, which provides a simple and universal picture of 
tunneling j20|. We then construct an effective action in 
D = 2 written in terms of composite fermion density and 
current, as well as the Chern-Simons gauge field. Sec- 
tion H ends by proving an important identity for this 
action which is used in the following part of the paper. 

In this paper we focus on the relatively simple "dirty 
composite fermion" case, corresponding to composite 
fermions scattered by the disorder, and described by fi- 
nite Ohmic conductivity. In Sec. Ill we consider the 
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problem with short-range interaction between compos- 
ite fermions. In the D — 2 action we integrate over the 
variables in the bulk and derive an effective D = 1 ac- 
tion that describes the dynamics in terms of the variables 
at the edge. This action is basically of chiral Luttinger 
form, with an extra "dissipative" term nonlocal in space 
and time, which takes into account the effects of charge 
relaxation in the bulk. The D = 2 — > D = 1 reduction 
for the problem with short-range CF interaction can be 
handled in an elementary way and leads to a simple al- 
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gebraic expression for the tunneling exponent in terms of 
Ohmic and Hall conductivities. 

Then, in Sec. |[y]we repeat the analysis for the problem 
with long-range Coulomb interaction. In this case the 
D = 2 — ► D = 1 reduction procedure involves solving 
a boundary value problem for dynamical screening near 
the edge. We consider three different models, describing 
the problems with unscreened Coulomb interaction and 
also taking into account screening due to image charges 
induced in the doped overgrown region. (This screening 
has the peculiarity that the screened interaction remains 
long ranged, because the image charges are located not 
above the 2DEG, but on the side of the 2DEG edge.) Two 
of these boundary value problems can be solved by ele- 
mentary methods using Fourier transform, and one leads 
to an integral equation of Wiener-Hopf type. In all three 
cases, we use the effective D = 1 action to compute the 
tunneling current, and derive an expression for the tun- 
neling exponent a. 

In the case of unscreened interaction the tunneling ex- 
ponent a turns out to be somewhat frequency dependent, 
having a contribution proportional to p xx In a;, which cor- 
responds to a slight deviation from the power law. How- 
ever, in the most realistic of the three models accounting 
for screening by the doped region, we find a nearly per- 
fect power law. Otherwise, the results for the three mod- 
els with long-range interaction, screened and unscreened, 
and for the short-range interaction model, give essentially 
the same dependence of the tunneling exponent on p xy , 
and thus all agree. The agreement of the results for dif- 
ferent kinds of interaction implies that they are robust. 

In the calculations described above, we characterize the 
system by a resistivity tensor that is independent of wave 
vector and frequency. In particular, this assumption im- 
plies that we are restricted to tunneling at voltages and 
temperatures small compared to the scattering rate of 
the composite fermions. At energies above the scatter- 
ing frequency, but below the Fermi energy, one is in a 
different regime (the "clean regime") where ballistic dy- 
namics should be used for the composite fermions. This 
regime may be of considerable practical interest because 
the samples used for the tunneling measurements are of 
very high mobility, and are presumably quite clean even 
near the edge. Even for electron energies below the CF 
scattering frequency, however, one should really check 
that contributions from wave vectors larger than the in- 
verse mean free path can safely be ignored. 

A proper treatment of the ballistic region requires the 
use of nonlocal electromagnetic response functions, and 
is considerably more difficult than the models discussed 
in the present paper. In the Appendix [A| below we inves- 
tigate a simplified model for the nonlocal conductivity, 
which serves to illustrate some of the salient features of 

I 



the problem. The simplified model is not adequate, how- 
ever, to answer unambiguosly the fundamental theoreti- 
cal question: whether low-energy degrees of freedom at 
short length scales can significantly alter the tunneling 
exponent at low electron energies. 

In order to better address this problem, we have also 
undertaken a numerical solution of the charge spread- 
ing problem with a proper representation of the nonlocal 
conductivity. Preliminary results suggest that the tun- 
neling exponents will not be changed by a large amount 
from the results obtained in the present paper [ p4| , but 
further work is necessary here. 

One should also recall that in the limit of very low 
temperatures and frequencies, in compressible states, one 
expects that there will be interaction corrections to the 
resistivity itself which depend logarithmically on energy 
0. Therefore, in principle, at sufficiently low energies, 
the renormalized value of p xx will become comparable 
to the value of p xy and our entire analysis may cease to 
be valid. However, the energy range where this would 
occur is too small to be of experimental interest in high- 
mobility samples where the bare value of p xx is small. 



B. The semiclassical phase method 

The tunneling exponent (|l|) was derived in Ref. |ic| l 
using a "semiclassicalphase" approach. Here we restate 
the derivation of Eq. (Q) emphasizing the connection with 
the effective action method being used in the main part 
of this article. 

One advantage of the semiclassical phase method em- 
ployed in Ref. is that it docs not require subtrac- 
tion of counterterms like S — iSf rcc used in the following 
sections. A suspicious reader may think of this subtrac- 
tion as an ad hoc procedure motivated only on physical 
grounds. Although we justify the counterter ms s ubtrac- 
tion carefully and rigorously below in Sec. II C, it will 



perhaps be useful for the reader to see the same result 
derived by an alternative method. 

It should be mentioned that the phase method, al- 
though more appealing intuitively, is more difficult in 
use, especially in problems with the boundary, like the 
edge tunneling problem. Because of that our use of it 
here is limited to the simplest case when the interaction 
is solely due to the Chern-Simons gauge field, and there 
is no long-range Coulomb interaction. The short-range 
interaction is assumed to be taken into account by the 
composite fermion transformation. 

We start with the tunneling electron Green's function 
in imaginary time. One can formally write it as an aver- 
age over the fluctuations of the gauge field: 



G rr i{t\2) = Z 1 / r D{a fl }G rr <(t 1 ,t 2 , a^)e" 



e//L 



t\2 — ti — ti 



(2) 
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r 



This exact expression emphasizes the order of integra- 
tion over fermionic fields and the gauge field a M . Here 
G Tr i{t\,t 2 , dfi) is the electron Green's function for a 
given configuration of the gauge field a M (r, t) . For eval- 
uating the tunneling current, we will need G rr <(ii2) for 
r = r'. 



The effective action £ 



mRef. §. 
order: 



is the RPA action derived 



Below we will only need C e f j up to quadratic 



a^(x)Vz};{x,x')a v {x')d 3 xd 3 x' (3) 



where the 

'Df/ l V\3') X ) " 

absence of long-range interaction in the RPA approxima- 
tion is given by 



correlator of gauge field fluctuations 
= (a fl (x) a v (x')) for the CF system in the 



Airp 



(4) 



Here K,^ v — (j^ju) is the free fermion current correlator 
(cf. Ref. § and Sec. below). 

We employ a semiclassical approximation for 
G rr (ti, t%, a M ). To motivate it, think of an injected 
electron which rapidly binds p flux quanta and turns 
into a composite fermion. The latter moves in the gauge 
field and picks up the phase 



Let us remind the reader that electron Green's func- 
tion in the composite fermion theory has an additional 

phase factor exp (j, ao{t')dt'^J introduced by Kim and 

Wen |23| which accounts for the gauge field of a solenoid 
inserted into the system upon the transformation of the 
tunneling electron into a composite fermion. This phase 
factor has been discussed in the context of the problem of 
tunneling into the bulk. By virtue of gauge invariance of 
electron Green's function under gauge transformations of 
the Chern-Simons field, one can eliminate the phase fac- 
tor using the Weil gauge a — 0. Because of that, seem- 
ingly different approaches to the bulk tunneling problem, 
some emphasizing the phase factor [ p3[ and others ignor- 
ing it 22 20(], are essentially equivalent. Below we are 
going to use the ao = gauge, which permits us to drop 
the solenoid phase factor from the start. 

Now, we substitute the Green's function in the phase 
approximation (g) into Eq.(|) and average over fluctua- 

This gives 



tions of a M using the action ( $ 



G„(t) = G<®(t)t 



-s 



where 



d 3 xd 3 x'f™(x)fi ee (x') V^(x,x'). 



(7) 



(8) 



a"(r, t)jjr(r, t)d 2 vdt, 



(5) 



where jjj rc0 (r, t) is the current describing spreading of free 
composite-fermion density. Semiclassically in a M (r, t), 
one writes 



G rr (ti,t 2 ,a fl )=e i ^ G^(t 12 ), 



(6) 



where G^(t 12) ~ i 12 1 is the composite-fermion Green's 
function in the absence of the slow gauge field. Note that 
fast fluctuations of a M are included in G^°\t) through 
rcnormalization of Fermi- liquid parameters. 
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i.e.. 



It is convenient to rewrite the exponent S hereafter called 
"action" as follows 



S = -- 



d 3 xf™{x)~a»{x) 



(9) 



where a^{x) —if T> liu (x, x')j^ cc (x')d 3 x' is the actual 
gauge field produced by the moving charge. The rep- 
resentation (g) follows directly from the ladder structure 
of the RPA response functions. 

From now on we adopt the ao = gauge, in which the 
relation between a and j takes the form 



a(r, t) = 47rp / z X j(r, t')dt' . 
J — 00 



(10) 



With this, the action S finally becomes 

s=£^(/«*)*j««*r) 



(ii) 



Note that we are working at T — 0, and the sum over 
Matsubara frequencies should actually be interpreted as 
j duj/2-K. From the form ( |Tl| ) we proceed to evaluate S. 



I 

The currents j^ eo (r) and j w (r) are found from the dif- 
fusion and continuity equations, 



j = -DVn, (w-VDV)n= J u (t) ; (12) 
= -l)(o)vn frec , (ui - V£> (0) V)n froc = J w (r), 



where J u (r) = e{e luJtl - e iult2 )8 {2) {r ~r G ). The diffusivity 
and resistivity tensors obey the Einstein relation 
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b-l 



a/3 — K Paf3, 



\ U >a/3 — K Pa/3 



(13) 



(0) i 
Pa/3 = p„p + 4?rp— e Q/3 



(14) 



where k is compressibility of free composite fermions. 
(Here "free" indicates the absence of long-range interac- 
tion, whereas the short-range interaction is assumed to 
be present and to give rise to the composite Fermi-liquid 
physics.) 

The resistivity tensors p and p^ are related by the 
composite fermion rule 

I 



We remark that, in our notation, the diagonal tensor 
elements of the imaginary time conductivities, resistivi- 
ties, and diffusivities have a sgnw dependence on w (see 
Sees. and [n| for details). Consequently, we may write 
D(u>) = ~D T (~uj) and n w (r) = -n_ w (r). 

Using these relations, one can simplify the expression 
for the action as follows: 



S = V 27Tp 

UJ 

,2 



d 2 r 



T — 

^ 2huj 

UJ 

Y.— 

UJ 

Ee 2 K 
2huZ 

UJ 

2 

E6 K 
2huZ 



(^i°lV Q ,n f - e (r)) e aP (B^V^r)) 
d 2 r (V Q n f - e (r)) (D^(p-p^)I)) (V^n w (r)) 
d 2 r (V Q n fre J(r)) (kD^ - kd) (V„n w (r)) 
d 2 r (n frc J(r)VI)Vn u (r) - ^(r)VI) (0) Vri f ™j 

(r) (-W + V£v) ^(r) - n^r) (-w + VflC'v) n frc J 



d 2 r 



_, free / 



E|£ / d 2 r(n f Z c (r)Ur) + J-UrM)) 
E / ^rJ-^r) K(r) - r£ cc (r)) 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



In the above equations, the tensors D and ZM°) are un- 
derstood to be always evaluated at frequency oj, not — oj. 
In going from Eq. (jlj) to Eq. ( |l8| ) we were able to discard 
the boundary term because the currents normal to the 
boundary are vanishing, as described below. The form 
( pl| ) will now be used for computing the action. 

The density n w (r) is found by solving the diffusion 
equation in the half plane y > 0, with the boundary con- 
dition j y — —DyydyU — Dy X d x n = at y = . In Fourier 
components n(x,y) = J2k elkXn k(y) this becomes 

(O 2 + q 2 )n k {y) = e(e^ - e^)6(y - y ), 

Dyyd y n k (y)y^ Q = -~iD yx kn k (0), (22) 



where q — (k 2 + u>/ D xx ) x ^ 2 ■ After solving this elemen- 
tary boundary value problem we take the limit yo — > 
and have 



D xx q + iDy X k 



(23) 



The expression for n free is similar, up to changing to 



y 

By inserting n and n h ' cc thus found into Eq.(|2l]) one 
obtains 



ui.k 



a iuJt 2 12 



&x%Q + ia- yx k a f) q (0) + 



(24) 



Note that this is precisely the expression for the action 
found in Ref. JhJ. Upon evaluating the integrals over k 
and u> it gives the result (0) in the limit a xx — > and a 



more general result (pi) for finite a x 



Note that the first term in Eq.(24), after integration 



over k and uj, is a smooth function of cr xy , whereas the 



G 



second term gives rise to a cusp in the tunneling expo- 
nent at a^y — 0, i.e, at v = 1/2. Indeed, the first and 
the second term of Eq.(24j) correspond to the first and 
the second term in Eq. (|l ) , respectively. This means that 
the plateau in the tunneling exponent for 1/3 < v < 1/2 
arises due to the second term. It is explicit in Eq.(|24|) 
that it is the second term that accounts for the free com- 
posite fermion dynamics, and so the cusp at v = 1/2 
should be understood as a signature of the composite 
fermion physics. 

Let us mention, that the expression (|24|) for the action 
can be rewritten as 



S 



2h 



1 



1 



u - VL>V 



(25) 



This formula can be taken as a hint that the problem of 
calculating the semiclassical action can be significantly 
simplified by a wise choice of an effective action and of 
a compensating counterterm. This is exactly what our 
strategy will be in the rest of the paper. 



as well as Chern-Simons fields. The tunneling consists of 
an instant process of adding one electron to the system 
and of its subsequent slow reaction. The second, coop- 
erative step involving Chern-Simons and Coulomb field 
relaxation controls the tunneling rate, while the first, 
single-particle step occurs instantly and contributes only 
to the prefactor in the tunneling current. Since for small 
bias the relaxation process occurs on a large scale, one 
may describe it using the semiclassical approach. How- 
ever, the fact that the tunneling particle obeys Fermi 
statistics is also important, and this will be included, fi- 
nally, in our analysis. 

In what follows we treat the system motion under the 
collective barrier semiclassically as classical Coulomb and 
Chern-Simons electrodynamics in imaginary time, find 
an instanton solution, and derive an expression for the 
tunneling rate in terms of instanton action. For that we 
generalize to the composite fermion system the semiclas- 
sical effective action theory introduced elsewhere p0[. 



B. Constructing the effective action 



II. EFFECTIVE ACTION IN D = 2 

A. Qualitative discussion 

Below we focus on the effect on tunneling arising due 
to relaxation of collective electrodynamical modes. Semi- 
classical theory can be used to describe it, assuming that 
the times and distances controlling the tunneling rate are 
large. 

The adequacy of the semiclassical approach can be un- 
derstood as follows. Tunneling in a strongly correlated 
system involves motion of a large number of electrons: 
While only one electron is actually transferred across the 
barrier, many other electrons are moving in a correlated 
fashion to accommodate the new electron. This collective 
effect becomes progressively more important as the bias 
decreases. At a small bias V, the single-particle barrier 
traversal time is much shorter than the relaxation time 
r ~ h/eV in the electron liquid. Therefore, while one 
electron is traversing the barrier other electrons essen- 
tially do not move. Thus instantly a large electrostatic 
potential is formed. The jump in electrostatic energy 
by an amount much bigger than the bias eV means that 
right after the one electron transfer we find the system in 
a classically forbidden state under a collective Coulomb 
barrier. In order to accomplish tunneling, the charge 
has to spread over a large area until the potential of the 
charge fluctuation is reduced below eV. If the conductiv- 
ity is small, the spreading over a large distance takes a 
long time, and thus the action estimated as the collective 
barrier height times the relaxation time r is much larger 
than h. 

This argument fully applies to a composite fermion sys- 
tem consisting of quasiparticles interacting via Coulomb 



The effective action can be written in terms of com- 
posite fermion charge and current densities n(r, t) and 
j(r, t), as well as the Chern-Simons gauge field a M . The 
total action is 



Stotal = Scf + Scs + $c 



Sb.c 



(26) 



In this section we motivate, define, and discuss different 
parts of the action ( p6| ) for our system. 

Below we focus on the case of diffusive CF transport 
taking place in the presence of disorder. Because the elec- 
trical conductivity is local in this case on scales larger 
than the mean free path, this problem is simpler than 
that of ballistic CF dynamics. 

The assumption underlying our analysis is that the 
main contribution to the action of the tunneling charge 
arises from large spatial and time scales, and thus local 
deviation from equilibrium is small. Therefore, one can 
expand the action in powers of charge and current den- 
sities, n(r,t) and j(r, t), and keep only the terms up to 
quadratic. 

The contribution ScF(n,j) is defined to correctly re- 
produce the equations of motion of composite fermions 
decoupled from the gauge field but interacting via 
the Coulomb potential. [To be more precise, since com- 
posite fermions describe interacting electrons in a mag- 
netic field, the short-range part of the Coulomb interac- 
tion is included in the definition of n and j of compos- 
ite fermions, so only the residual long-range part of the 
Coulomb interaction enters the action ScF(n,j).] We 
consider ScF(n,j) of a quadratic form constructed using 
CF response functions. One can see that the requirement 
of matching the CF equations of motion is not entirely 
sufficient to determine the action, e.g., because it leaves 
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freedom of rescaling the whole action or even its differ- 
ent parts corresponding to different normal modes of the 
problem. 

The exact form of the action can be determined in the 
following way pc| ] . The action used to study the dynam- 
ics in imaginary time is precisely the one that appears 
in the quantum partition function. The latter action ex- 
panded up to quadratic terms in the charge and current 
density must yield the correct Nyquist spectrum of equi- 
librium current fluctuations: 



7 af3\ 



'D/3p'q a >qp> 



(27) 



Here g(r) = j(r) + D\7n(r) is the so-called external cur- 



J 



CF 



-Y 

2 ^ 



rent. In this article we are interested in the hydrody- 
namical regime of small frequency to and momentum q, 
in which case the conductivity and diffusivity tensors a a p 
and D a p satisfy the Einstein relation a — c 2 kqD, where 
ko = dn/d^i = to* /2irh 2 is the free CF compressibility. 
Generally, both a and D are functions of to and q. 

Below we assume isotropic conductivity tensor charac- 
terized by a xx and a xy . Also, to make expressions less 
heavy, we often use the units h = e = 1 in intermediate 
steps of calculation, and recover h and e in final results. 

The requirement of matching equilibrium current fluc- 
tuations is essentially equivalent to the fluctuation- 
dissipation theorem. The action in imaginary time reads: 



d 2 rd 2 r' 



g 1 u (r)K a p(u,, r, r')g^(r') + U(r - r')n_„(rK(r') 



(28) 



where U(r— r') is the electron-electron interaction, and the kernel K(uj, r, r') is related to the current-current correlator 



S°) (, ,\ n(°) 



(29) 



given by Eq.(§7|). Here ct (0) (w) and are func- 

tions of the Matsubara frequency ui obtained from the 
real frequency functions by the usual analytic continua- 
tion. The superscript (0) here and below indicates that 
the response functions cr' ' and correspond to the 
free CF theory, in the absence of coupling to the Chern- 
Simons field and interaction U(ri — r2). 

It is appropriate to recall here the general properties 
of the Matsubara conductivity a a p (uj) . By the symme- 
try of kinetic coefficients, the dielectric function is an 
even function of the Matsubara frequency: e Q( a(iw) = 
£/3a(— iuj) (see Ref. Q). Relating it to conductivity by 
ea/3(w) = S a f3 + 4-7T a a p(oj)/iuj, one obtains that the lon- 
gitudinal (Ohmic) conductivity is an odd function of lo, 
while the Hall part is an even function of uj. This means 
that the constant conductivity case actually corresponds 
to a discontinuity in <t xx (lu) at uj — 0: 



,(iuj) 



a xx sgnw 



(30) 



whereas a xy (iuj) = a xy has no discontinuity. The 
same applies to the components of the diffusivity tensor 
D a p(uj). 

The coupling of composite fermion charge and current 
to the statistical gauge field a M (r, t) is described by the 
Chern-Simons action in a standard way : 

S C s = i J dt J d 2 r(na + j • a + — e^a^ax 

(31) 



Here p is an even integer corresponding to the number of 
flux quanta in the construction of composite fermions. 

The charge and current densities entering Eqs.(p8|) and 
( |3~i| ) are not independent. They may satisfy a continuity 
equation. For tunneling problem we employ 



Vj = J(r,<) 



(32) 



where the source J(r, t) = e5(r — ro)[S(t — t\) — S(t — £2)] 
describes adding a composite fermion at the time t\ at 
the point ro and subsequently removing it at the time ti 
at the same point. To handle this constraint, one has to 
put in the action (p6|) the term 



Scant = i J (n(r,t) + Vj(r,t) - J(r,t)) $(r,t) d A rdt. 

(33) 

with the Lagrange multiplier function $(r,t). 

Finally, to complete the action, one has to ensure 
proper boundary conditions. We choose the coordinates 
so that the 2DEG occupies the half plane y > 0, so that 
the half plane edge coincides with the a;-axis. The bound- 
ary conditions at the edge arise from the requirement that 
normal current at the edge vanishes: 



j y (x,y = 0,t) = 



(34) 



The corresponding part of the action is constructed by 
using another Lagrange multiplier: 



5b. c. = i dx dt j v (x,y = 0,t)cf)(x,t) 



(35) 
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Besides ensuring proper boundary conditions at y = 0, 
the term (|35|) is needed to make the total action gauge 
invariant with respect to gauge transformations of the 
Chern-Simons field a u . 

above, we do not need to in- 



As remarked in Sec. 



IB 



elude in the effective action a term expressing the effect of 
the solenoid that appears in the system upon the trans- 
formation of the electron into a composite fermion. Since 
we will work in the ao = gauge, the "string" phase fac- 
tor exp (i J t * 2 ao(t')dt'J of Kim and Wen Q is absent. 

As a validity check of the action (^6|) let us derive the 
dynamical equations. They are obtained by taking the 
variation of the action (|2^) with respect to all variables 
excluding the Lagrange multiplier <E>(r,i). The resulting 
equations are of the standard form: 



2np 
1 



P< 0) j = E CS - 

E c \ s =s<*< 3 f 



VUn 



Bcs =n + J 



(36) 
(37) 

(38) 
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where ~Ecs = Va + a and Bcs — V x a are Chern- 
Simons electric and magnetic fields. The effective inter- 
action U is defined as 



1 



U(r - r') = U(t - r') + — S(r - r') 



Kq 



(39) 



where U(r — r') is the electron-electron interaction and 
Ko = ra^jlixh 1 is the compressibility of free composite 
fermions. Both U and p(°) in Eq.(|36|) in general act as 
nonlocal operators. The boundary condition j y — 0, ac- 
cording to Eq. (|37|) , requires that the tangential Chern- 
Simons electric field vanishes at the boundary: d x = 0. 

Also, it is straightforward to check that eliminating the 
Chern-Simons field leads to Ohm's law with a corrected 
resistivity tensor: 



pj = — V(Z7n), where p — p^ ' + 



ph 



-1 

1 



(40) 



is the measured resistivity tensor. Note that Chern- 
Simons interaction changes p xy , while p xx remains intact. 



C. The fundamental identity 

The nonlocal current-current term in Eq.(|28|) makes a 
calculation for the problem in the half plane y > long 
and not too transparent. To circumvent this algebraic 
difficulty, we derive an identity for the action (Eq) that 
allows us to replace it by an equivalent action with a local 
current-current term. 

To that end, we introduce another CF action: 



^loc 
'CF 



2 ^ 



d 2 vd 2 v' 



~3-Mp%{^ r')£(r>) + U(r r')n_ w (r)n w (r') 



(41) 



where uj is Matsubara frequency. Here p^l(v, r') = 
p^gS(r — r') is the resistivity tensor and U is defined 
by Eq.©. 

The relation between the actions ( pq ) and (|4^) is pro- 
vided by the following fundamental identity: 



ScF{n,j) = Scp(n,j) 



$C F i. n frcc , Jfrec) , 



(42) 



where n(r,t) and j(jc,t) are arbitrary functions satisfy- 
ing the continuity equation ( |32"| ) and the boundary con- 
dition (|3^), whereas nf rco (r, t) and jf lee (r, t) correspond 
to the saddle point of the action describing noninteract- 
ing composite fermions decoupled from the gauge field. 
Thus the functions 7if ree and jf Tee can be found by solving 
Eqs. ©-© with U(r-r') = ^^(r-r') and no E cs 
and Bcs- Supplemented with the continuity equation 
that is present in the effective action ( |26| ) as a constraint, 
the equations for fif rco and jf rco take the form: 



jfree(r,u;) = -5 (0) Vn frcc (r, uj); 
cjn fre c(r, uj) + Vj froo (r, uj) = J (r, uj) 

The boundary condition for the system (^~; 
sence of normal current jfree 

at y = 0. 



(43) 

t) is the ab- 



The result ( f42| ) is formulated and established below 
for local resistivity, because in this case the proof is more 
straightforward. It is possible, however, to generalize it 
to the case of nonlocal resistivity p^(r,r'). This re- 
quires more general arguments which will be discussed 
at the end of this section. 

To prove the identity fl42] ) , we write the expression (|2 
for the kernel K~ x using gradients: 



K- 



.a/3 



(44) 



where the operator convention is that ^ Q acts to the 
right, whereas ^ Q acts to the left. It is useful to intro- 
duce the distinction between V and ^ and to keep track 
of it later, so that we are able to invert the kernel K~ x 
and to evaluate the expression in the first term of the 
action ( p8| ) before doing the integral over the half plane. 
In this way we can properly handle boundary terms. 

Inverting Eq.([4~4|) and using the Einstein relation be- 
tween D a p and a a p together with the relation between 
conductivity <j a p and resistivity p a/ 3, one obtains 
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K a = — ^ - \? r (45) Consider the first term in the action (p8|): 

I 



S-uK a ^Z = -S-Md - (V ■ g-u,) m (V • g u ) (46) 

w k w(w + VDV) 

Below we perform some manipulations with the expression (fiff), refraining from integrating over r until the very end, 
because of the above-mentioned need to be careful with gradients and boundary terms. 
Now we substitute 

gS=j° +D a P(Lu)Vpn (47) 
in the first term of the right-hand side (RHS) of Eq.(Q), and find 

-S-Md = -3 a p a pf + — (j a V a n + (V a n)j a ) + — n (VQV) n. (48) 

LU LU KqLU KqUJ \ / 

To transform the second term of the RHS of Eq. (^6|) , we substitute 

V • g = V ■ j + ujn - (uj + *VD~^)n (49) 

and obtain 

- (V • g_ w ) — ^=j=- (V ■ g w ) = -J V^n - 7 + —{Jn + nJ) - — n L + VBV) n, (50) 

where J = V • j + am. 

Finally, we add the expressions ( jig ) and (|50|), and combine the last term in Eq.(^) together with the second and 
third terms of Eq. (|5(|) . After doing this we find the resulting expression 

Vpa/?/ + — n 2 - J ^ + — V Q (j Q n + (51) 



Upon integrating this expression over r and multiplying 
by 1/2, the first two terms give corresponding terms of 
the action (fill) , the third term gives <Sj?£. (7if ree , jfree) ap- 
pearing in (|42j), and the last term vanishes due to the 
boundary condition (j^), thus proving the identity d42|) 
Having given a formal proof of the identity fl42|), let 
us now point out the relation of Eq. (fl2|) to the struc- 
ture of RPA diagrams in the perturbation theory for 
Green's functions in the presence of disorder. To sim- 
plify the discussion, let us ignore the CS gauge field, 
and consider the problem of electrons coupled only by 



Coulomb interaction. In this case, the RPA self-energy 
£ can be represented graphically, as shown in Fig. |]. 
In the D = 2 problem the bare unscreened interac- 
tion, represented in the figure by a thin broken line, is 
£/(k) = 27re 2 /e|k|. The diffusive polarization operator 
is n(k, to) = kqD\l 2 /(lu + £>k 2 ), and the diffusive vertex 
part is 1/(uj + Dk 2 ). One can verify, by performing a 
resummation, that the dynamically screened interaction, 
shown in Fig. [l] by a thick black line, can be represented 
as follows: 



1 



C/(k) 



(cu + Dk 2 ) 2 l + U(k,Lu)U{k) 



Dk 2 



+ l + K^tk) 



L>k 2 



(52) 
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as the difference between the propagator of an auxiliary 
interaction and the diffusive vertex part, multiplied by 
lo^ 1 . These two contributions are shown in Fig. [l]by the 
wavy red and wavy green lines, respectively. 

The self-energy diagram in Fig. ^ corresponds to in- 
teraction via a dynamically screened Coulomb potential, 
i.e., to a shakeup of a dissipative plasmon. This effect 
is described by the hydrodynamical effective action in- 
troduced above in Sec. Q, and so it is to be expected 

I 



that the expression in the RHS of Eq.(p2[) corresponds 
directly to the difference <S loc - S£ c c c in Eq.©. 

On can rewrite the formula (|5|) in a quite general op- 
erator form, generalizing it for any interaction U , polar- 
ization operator II(w), and vertex part V(u>) 7 satisfying 
the Ward identity II (u>) = Kq(1 — ljV(u>)). For that, one 
represents the vertex part in the form V(v) = (X + uj)^ 1 , 
and writes: 



v(u>) (l + c/nH 



1 UV(u>) = - 



1 



^ V X + UJ (1 + kqU)- 1 X + lu 



(53) 



The formula ((53|) is proven straightforwardly, by expand- 
ing the fractions in operator geometric series, and subse- 
quent resummation. 






+ + o + - = 



V 



A 



FIG. 1. Resummation of the RPA diagram series for 
self-energy for diffusive electrons. The black wavy and broken 
lines represent dynamically screened interaction and bare un- 
screened Coulomb interaction, respectively. The bubbles and 
triangles represent diffusive polarization operators and vertex 
parts. The red and green wavy lines on the RHS are defined 
in Eq. E| 



One can view the formulas ( |52| ) and (|53J) as a motiva- 
tion for the identity (fl2|). More importantly, the relation 
to RPA diagrams, explicit in Eq.(p2[) and (|5^, demon- 
strates the general character of the identity Q42[) , which 
is not evident from the way it is justified above. Compar- 
ing to Eqs.(52) and ( p3| ) makes it clear that the identity 
( [4 4 ) is robust under changes in the geometry of the sys- 
tem, alterations of the boundary conditions, and addition 
of more complicated interactions such as the CS gauge 
fields. 

I 



I 

The analog of Eq.(^) and Eq.([53"|), and thus of the 
identity ((42|), holds even for ballistic Fermi-liquid dynam- 
ics. In this case, according to the microscopic theory of 
Fermi liquid, II = K kv/(kv — to) and V = l/(ui — kv), 
and the operators act on the particle-hole distributions 
on the Fermi surface. For a Fermi liquid, the formula 
(p3|) holds with A = kv. 



III. THE D = 1 ACTION FOR SHORT-RANGE 
INTERACTION 

A. Integrating out variables in the bulk 

In this section we consider the simplest model of short- 
range interaction, U(r — r') = US(r — r'), and diffusive 
CF transport described by p a p{r, r') = p Q /3<5(r — r'). 

We shall start with the action Stotai given by Eq.(|2^) 
in the half plane and derive an effective D — 1 problem 
by integrating out the dynamics in the bulk, and keeping 
only the variables at the edge. Since the action (26) is 
quadratic, the integration can easily be performed by the 
saddle point method. 

From now on we replace the CF action ( p8| ) by the ac- 
tion ( f4l| ) with a local current-current term. The virtue 
of doing this is that the action ([0]) is much easier to 
handle, whereas the identity (|4|) allows us to go back to 
the physically meaningful action (|2^) at the very end. 

First, it is convenient to integrate out the Chern- 
Simons gauge field a M , both in the bulk and at the edge. 
We do it by fixing the gauge ao — 0. Upon integra- 

into 



Pxx — 

with 



Pal turnS 
- p$j + ph/e 2 

p^x ■ The action acquires the form iS to t = S — Sf Te , 



tion over a M the CF resistivity tensor 
the electron resistivity tensor (fh 



'xy 



s 



^-n_ w (r)n w (r) j + S cont + S hj 



(54) 
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Then we integrate out n and j in the bulk, keeping fixed the normal current j y at the edge. The result is 



S = J2j d2r (l wa a/3 (o;)V Q $_ w (r)V^$ a; (r) + ^^(r^r) + i$(r, t) J(r, t) 

+ i J dxdt [&(x, y = 0,t) — 4>{x,t)] j y (x, y = 0,t). (55) 

Here <t(cj) = p _1 (cj) is the electron conductivity tensor. The frequency dependence of &(uj) is the same as that of 
p{u)): (t xx (uj) = cr^sgncj, a xy (u>) = <j xyi etc. 

The next step is to integrate over j y (y = 0), which gives <&(x, y = 0, t) = 4>(x, t). Hence, the action is 



S = J2j ^ 2 ^Q^a/3(w)V Q $_ w (r)V^^(r) + ^$_ w (r)$ w (r)^ +i J dxdt® 



(x,y = 0,t)J(x,t). (56) 



In handling the source term J we assume that the point 
fo = (xo,yo) at which charge is injected is very close 
to the boundary, i.e., yo — > 0, and thus the source in 
Eq.(^6|) can be effectively placed at the edge: J(x,t) = 
ed(x-x )[S{t-ti)-d(t-h)}. 

Finally, we integrate out the bulk value of <E>(r,i). 
From Eq.(p6|) the equation for $ at y > is 

a xx (Lu)WHUr) + |* w (r) = 0. (57) 

It is convenient to use the Fourier transform of $ w (r) 
with respect to variable x only: 

= X><-,*(tf)e <fcB . (58) 
fc 



Note that Fourier transform in y is not suitable because 
we are dealing with the boundary value problem in the 
y > domain. 

Then the solution to the equation ( |57j ) is straightfor- 
ward: 

* u ,k(v) = *»AV = 0)e- q(u ' k)v ; (59) 
q 2 (u,k) = k 2 + 

U <T XX 

After substituting Eq.(]59]) into (|56|), one obtains a D = 1 
action: 



5 = - (cr xa: |w|g(a;, fc) + ia xy u>k) fe + J(-w, -fc) . 



(60) 



where we put Eq. ( p0[ ) in the Luttinger liquid theory form 
in terms of the boundary field </>(x, i) = ^(a;, y = 0, t) in- 
troduced above as a Lagrange multiplier. 

This effective action represents a generalization of the 
chiral Luttinger theory of edge modes to the compressible 
problem with finite a xx . Because of the relation between 
q and u>, the dissipative term in the action (|6^) is non- 
local in the time representation. In the incompressible 
limit a xx — > 0, we recover the standard chiral Luttinger 
action: 



B. Instanton action 



S = — I d x 4> d t 4> dx dt + i 



J(x,t)4>dxdt (61) 



In the above derivation we ignored effects of the boundary 
compressibility. Taken into account, these effects lead to 
an additional term of the form j (d t (f>)' 2 'dxdt which does 
not affect the long-time dynamics and drops from the 
final answer for the instanton action derived below. 



The source term in the action (60) describes coupling 
of the tunneling charge to the field <fi(x,t). Thus, the 
electron creation operator can be written as ip + (x,t) = 
tp^ F (x,t)e le ^ x,t \ where ipg F (x,t) is the operator of a 
composite fermion, and e is the electron charge. Let 
us point out the resemblance of the exponential e le $( x -t) 
to the standard one-dimensional Luttinger liquid expres- 
sion. 

Tunneling is related to the electron Green's function. 
To find the tunneling rate, we evaluate the equal point 
Green's function G(t) = (ip(0, ti)tp + (0, t 2 )} T =t 1 -t2 of an 
electron. Using the above relation of tp and ipCFi we 
write the electron Green's function in terms of the CF 
operators and then make a factorization approximation: 
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(V(0,*i)^ + (0,i 2 )) 



{i/>CF(0,h)il>£ p (p,t2))(exp (i 



J(x, t)<p(x, t)dxdtj 



(62) 



where the first and the second averages on the right-hand 
side are taken over the fermionic ground state and over 
fluctuations of the electric and CS gauge fields, respec- 
tively. This approximation holds because the dynamics 
of the injected quasiparticle and of the collective charge 
relaxation mode are decoupled in space and time. The 

I 
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CF quasiparticles and edge magnetoplasmons differ both 
in the rate of penetration into the 2DEG bulk and in the 
velo city of motion along the edge (cf. the discussion in 

Thus the imaginary time Green's function can be writ- 
ten as 



G(r) = G cf (t) (exp (i J J(x,t)<f>(x,t)dxdt) ) = G CF (r) e -[SM-s fcco (r)] _ 



(63) 



where Gcf{t) is the Green's function of a free com- 
posite fermion injected and later removed at a point 
of the boundary. In the last term of Eq.(|65|) we used 
the identity ([42] ) relating the average over cf>(x, t) in 
(exp (i J J{x,t)(j){x,t)dxdt)) to the action (p0[). 

According to the CS Fermi liquid theory, in the effec- 
tive composite fermion mass approximation, Gcf(t) — 
1/t. This essentially free fermion result holds even 
though the gauge field fluctuations give rise to infrared- 
divergent logarithmic corrections to the effective 
mass m*, because these corrections are canceled by cor- 
rections to the residue Z of the Green's function. 

The tunneling current is obtained from G(r) in a stan- 
dard way. One has to continue G(r) from imaginary to 
real time, and to do the integral over time: 



I(V) - Im / G(t) 
Jo 



JeVt 



-dt. 



(64) 



Now, we evaluate (exp (i J J{x,t)4>{x,t)dxdi\) using 
the local action (|60|). By a Gaussian integration, the 
result is e~ s , where 



|JHI 5 



2 ^ \u}\{a xx q + ia xy k sgnw) 



(65) 



The substitution k = fcosinh2a; with fco — (M/er^t/) 1 / 2 
simplifies integration over k: 



\J{lo)\ 2 f(e x + l)dx 



ere- 1 



(66) 



J 



I 

where a — a xx + iu xy . The integral ( |66| ) is taken in 
the domain — x max < x < a; max , and gives an ultravi- 
olet logarithmically divergent answer which we cut at 



s = 



M 



Px X , 4:k max (T xx U 
8tt 2 n 



4tt 2 



PxyVH 



(67) 



Note that this expression does not vanish even in the 
absence of interaction with the Chern-Simons Jreld and 
electron-electron interaction, when p = and U = Kq 1 . 
This indicates that part of the answer represents the con- 
tribution of noninteracting composite fermions and must 
be subtracted off. This subtraction happens automati- 
cally because of the identity (f42|), which confirms that 
the correct action is indeed S — Sf lcc . 

One can see that the counterterm <Sf rcc is indeed re- 
lated to the effect of free composite fermions. The phys- 
ical origin of the ultraviolet divergence at A; max is that 
for free fermions the relaxation is fast and involves large 
momenta k ~ kp. On the other hand, the contribution 
resulting from the interaction should not diverge at large 
momenta. 

To find S — <Sf roc , we subtract from Eq.(|6"7j) the same 
expression with p = and U = n^ 1 . Integrating the dif- 
ference over uj, we get S — <Sf rcc = (a — 1) lni/io, where to 
is a microscopic time of the order of the scattering time, 
and a is given by 



2e 
■nh 



2 r 



e 2 p xx 
irh 



In 



(1 + K U)a xx /cr 



(68) 
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where Oh = tan 1 p xy j p xx is the Hall angle, U is the 
short-range interaction, and k$ — m*/2iT?i 2 is the free 
CF compressibility. The behavior of a as a function of 



D 



a 
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FIG. 2. The tunneling exponent for the model (|3£|) with 
short-range interaction U(r — r') = US(r — r'), where 
k = Ko = m,/2nh 2 . 

To verify that a is the tunneling current exponent, we 
write the electron Green's function as Eq.(63), where the 
free composite fermion Green's function is Gcf(^) ~ t -1 - 
Therefore, the Green's function is G(t) ~ t~ a . One can 
compute the tunneling current from Eq.(j64|), and obtain 
the power law I(V) ~ V a . The expression ( |68|) shows 
that the shake-up effects suppres tunneling in a uniform 
fashion for the filling factors v both on and off the quan- 
tum Hall plateaus. The I — V curve is given by a power 
law with the exponent depending smoothly on the fill- 
ing factor, via the resistivities p xx and p xy , and effective 
interaction kqU. 

One can compare this result with the chiral Luttinger 
liquid theories of tunneling into the edge of an incom- 
pressible QH state. For that, one has to consider the 
limit of a large Hall angle: 9h = #2 = tt/2. In this 
case p xx — > and the exponent ([38]) acquires the form 
(|l|) corresponding to a staircase with plateaus in the in- 
tervals 1/3 < v < 1/2, 1/5 < v < 1/4, etc., interpo- 
lated by straight lines with the slope 2. At the rational 
filling fractions v = n/(pn ± 1) we recover the results 
of the Luttinger liquid theories. To see this, substitute 
Pxy — (p + l/n)h/e 2 , p x °y = h/ne 2 in the expression 
(0), and get / - yi+\p+i/n\~i/\n\ > which agreeg with thc 

universal tunneling exponents predicted by Wen and by 
Kane, Fisher, and Polchinski for Jain filling factors with 
positive and negative n. 

It is interesting that the tunneling exponent (Q) has 
cusplike singularities near the compressible rational u's 



with even denominator, v = 1/2, 1/4, etc. The origin 
of this effect is a qualitative change in the structure of 
the edge modes near these filling factors. In particular, 
let us discuss the vicinity of v — 1/2, where the quantum 
Hall state can be described as a Fermi liquid of composite 
fermions carrying p — 2 flux quanta each, and exposed 
to "residual" magnetic field SB — (2-^ _1 ). At v < 1/2 
the residual field direction coincides with the total field, 
and all edge modes propagate in the same direction. On 
the other hand, at v > 1/2, the structure of the edge is 
qualitatively different, consisting of modes going in op- 
posite directions. This effect makes v = 1/2 a singular- 
density from the point of view of the tunneling exponent. 

The singularities at v = 1/p are smeared in the pres- 
ence of scattering by disorder, i.e., at finite p xx . Inter- 
estingly, the deviation from the staircase described by 
the expression (|l|) due to effects of finite p xx can be ei- 
ther positive or negative, depending on the interaction 
strength k,qU (see Fig. |^). In the absence of interaction, 
at U = 0, the tunneling exponent a < a(p xx = 0). On 
the other hand, at large interaction, a > a(p xx = 0). 

It is instructive to compare the results and (|l|) 
with the exponent a = \jv found using hydrodynamical 
approaches [^- ^|Jl^ , [l7|l in which the edge dynamics is 
modeled as a charged fluid, without any additional in- 
ner quasiparticle degrees of freedom. Our expressions 
( |68| ) and (|l]) have the form of a difference of two contri- 
butions, the first of which is essentially \jv with small 



corrections due to finite 



Thc second contribution is 



expressed in terms of the response functions of free com- 
posite fermions, and it is this term that leads to nonan- 
alyticity and plateaus in a{y). According to the identity 
(|42|), these contributions arise from the local action 5i oc 
and the counterterm S^. e , respectively. It is easy to see 
that there is a direct correspondence between our action 
Sioc and the hydrodynamical actions |ll| 0,|l6],|l7| . In 
our approach, the role of the counterterm <S^ C C is to en- 
sure that the Green's function of free composite fermions 
agrees with Fermi statistics. From that point of view, 
the plateau- like structure in a(y) is a manifestation of 
the role of composite fermions as underlying quasiparti- 
cles of the QH state. 



IV. MODELS WITH A LONG-RANGE 
INTERACTION 

A. The action for the edge mode 

We assumed above that the interaction has a short 
range. Due to the long-range character of the Coulomb 
interaction, electromagnetic modes in a real system are 



very different from those considered in Sec. [II. Hence 
the effect of shakeup of these modes on tunneling is also 
somewhat different. In this section we extend the method 
outlined above to the problem with Coulomb interaction, 
and consider several situations describing screening of the 
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interaction in the overgrown cleaved edg e sy stem, as well the effective action for the edge outlined in Sec. flf can 
as the unscreened Coulomb interaction |^5| . be followed without any change up to Eq. (|55|), which in 

For the long-range interaction, the method of deriving this case takes the form 



J 



S = \ / d 2 r [^(wJV^WVjJ^r) + i$(r, t) J(r, t)] 
1 u J y>° 

+ I fd'r<Pr'(J*$- u (T)U- 1 (T,r')$ u (T')+i [ dx dt [$(x, y = 0, t) - <f>{x, t)] X (x, t), 



(69) 



where U 1 (r, r') is the inverse of the interaction kernel, 
and the notation 

X(x,t) = jy(x,y,t)\ y=0 (70) 

is introduced. It will be convenient now, instead of inte- 
grating over j y (x, y = 0, i) as we did above, to keep it as 
a dynamical field. 



Let us note that in the interaction term in Eq.(69) the 
integral over r and r' goes over the whole plane, not just 
over the half plane y > as in Sec. [II. The reason is 



simple to understand by writing the relation between $ 
and n: 



w$(r) 



U(r, r')n(r')(2V , 



(71) 



y'>o 



and observing that for long-range U the field $(r) ^ 
for both y > and y < 0. 



I 

To proceed with deriving the effective D = 1 action, 
we decompose the conductivity tensor into the diagonal 
and off diagonal parts, o a p((jj) = a xx agnu> SaB+ a xy e a p- 
The off diagonal conductivity term in Eq.(p9) is a full 
derivative, because 



e^V a $_ u (r)V^ u (r) = V Q [e^-^V^'Mr)] 



(72) 



As a consequence, this term is converted into the bound- 
ary term expressed in terms of $> y= o(x,t) 
the total action can be written as 



"total — 



2D 



■Si 



D ! 



where 



(x,t), and 
(73) 



J 



Sid =i I ( -(r X yd x (j)*{x, t)dt<f>(x,t) + (j){x,t)J(x,t) + [§ y=0 {x. t) - (f>(x,t)] x(x,t) ) dtdx 



(74) 



and 



S2D = \Y,(J y Q \u\a X xV a $- u (r)V a $ u (r)d 2 r + J J w 2 $_ w (r)f7- 1 (r,r')$ w (r')d 2 rd 2 r' 



(75) 



We included the source term J in Sid by placing it at 
the boundary y = and accordingly added the term 
i(f)(x,t)J(x,t) to Eq. ([74]) , simultaneously removing the 
term i$ u (r) J_ w (r) from Eq.(|75|), 

Now, one can integrate over the field $ w (r). This 
amounts to taking the saddle point of Stotai, i- e -' to solv- 
ing the problem 

-Mcr^V 2 ^!-) + wn(r) = iXu>{x)5(y) (76) 
w$(r) = / U(r,r')n(r')d 2 r' 

Jy>0 

in the domain y > with the boundary condition 
<9 y < I > !/ =o = which corresponds to the absence of cur- 
rent normal to the edge. This problem describes the re- 
sponse of the charges in the conducting half plane to the 



external charge source Xu{x)S(y). The solution of this 
problem taken at the boundary y = can be written 
as some linear operator applied to the source Xu) (x). In 
terms of Fourier components one has 

%=o(k,Lo) = Q- 1 (k^)I^L i (77) 

\W\(Txx 

which defines the function Q(k,uj) playing the key role 
in what follows. Interestingly, there is no dependence in 
the problem ( |7^ ) on a xy whatsoever, because the corre- 
sponding part of the action is a boundary term, and thus 
it belongs to the boundary action (|74"|). 

We postpone the discussion of the problem ( |76| ) 
and proceed with deriving the effective D = 1 ac- 
tion. The integration over <& w (r) simply adds the 
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term 5 X)fc ^ M 1( 3 1 {k,uu)x-k-uXk,u to the action Finally, we integrate over the field x( x > an d obtain 



Sid given by Eq.u[4 



the total action in terms of the boundary field (j>(x,t): 



S = ^ 9 l^xx M Q(w, fc) + i(T X yU}k]<fr-k,-u><t>k,u + <f>-k,-uJk 



(78) 



This action, in which the function Q(w, fc) has to be found 
by solving the problem ( [76]), r epresents the analog of the 
action (J60|) derived in Sec7|l| for short-range interaction. 

Using this action for calculating the Green's function 
goes in a complete parallel with section III. The re- 
sulting Green's function is G(r) = e _s G^(r), where 

g[?p(t) = r _1 is the free CF Green's function. The sad- 
dle point action S, by virtue of the identity ([42|), can be 
written as S = S loc — Sj° e c e , where <S loc and <Sf° e c e are found 

I 



I 

by taking an appropriate saddle point of Eq.([78|). The 
result is conveniently expressed in terms of a "spectral 
weight" A(to): 



1 



G(t) = — exp 

T 

J(w) = e(l - e- WT ). 
Here A(u>) is defined as 



\J(u)\ 2 A(u) — ■ ■ 
11 47m|cj|/ 



(79) 



a xx Q(k,Lu)+ia xy k a x °jQ( )(k,u) + 




(80) 



where Q(fc, w) is defined by Eq.(77), and Q^(k, lo) is de- 
termined from Eq.(f76|) for U(r, r') = Kq 1 S(t — r'), which 
corresponds to noninteracting composite fermions. While 
deriving fl80|), we replaced a xy ujk by a xy \ui\k in the ac- 
tion (78), which does not change the integral in Eq.(|80|) 
because a sign change of lu can be accommodated by a 
sign change of k. 

The relation between the tunneling exponent a and 
the spectral weight A(u>) is most simple when A does 
not depend on w, as in the case of short-range interac- 
tion discussed in Sec. III. In this case, simply a = A+ 1. 
A frequency dependent A(lo) can be interpreted as an 
energy dependent tunneling exponent 



a(u) = A(u) + 1. 



(81) 



This interpretation is meaningful only if the in- 
dependence of A is sufficiently weak. This will turn out 
to be precisely the case below, for the problem of long- 

I 



model V —V' 
model V + V' 
model V 



t/(r,r') = 
£7(r,r') = 
U(r - r') 



range Coulomb interaction, in which A{u>) varies with uj 
not faster than logarithmically. 

In what follows we consider the problem (|76|), find 
Q(k,oj), and evaluate the spectral weight ([so|). 



B. Solving for Q(u),k) 

The problem ( [76] ) that has to be considered in order to 
find Q(u},k) involves a long-range kernel U(r,r') and, in 
general, requires solving an integral equation. This equa- 
tion is defined in the half plane y > 0, and thus cannot 
be treated by simple tools. Generally speaking, one has 
to treat it by the Wiener-Hopf method. 

However, there are special cases corresponding to in- 
teraction screened by a mirror image in the region y < 
that can be handled by the Fourier transform. Below we 
consider three models: 



\ 8(r — r ) ; 

e|r — r'| e|r — r"| ko 

+ -6(r-v>) ; 



e|r — r'| e|r — r"| kq 

.2 1 



er-r' k 



(82) 
(83) 
(84) 
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Here the point r" is a mirror image of r' with respect to 
the edge y — 0: r' = [x 1 ', y'), r" = (x 1 , —y'). 

We start with the model V — V because it is sim- 
pler, and also because it directly corresponds to the over- 
grown cleaved edge system where screening of the type 
( p2| ) occurs due to the charges induced in the doped re- 
gion. One can transform the problem ( |7^ ) in the half 
plane y > to a problem in the full plane by ex- 
tending the functions n, and X to the negative half 
plane y < with a sign change: <fr(x, — y) = —&(x, y), 
n(x,—y) = —n(x,y). Similarly, the source X m ( K§ ) 
must be extended so that Xu (x, ~ v) = — Xu( x :V)- in 
that, the source Xu{ x iV) lSi assumed to be located not 
right at the line y = but somewhat away from it, so 
that the dependence of x m (HI ) below on y is given 
by Xu>(3, y) = Xu(x) [S(y - y ) - S(y + y )} with a small 
ya > 0. The limit yo — > will be taken at the end. 

Upon extending the problem to the whole plane the 
interactio n (|S2"|) has to be replaced by the unscreened in- 
teraction (|S4|). Then the problem ([76]) takes the form 

(-\uj\g xx V 2 + uj 2 U^ $ w (r) = A5\y) + i Xu ( x , y), 

(85) 

where U 1 denotes the inverse of the operator with the 
kernel @. 

The term AS'(y) is inserted because the function 
$^(r), extended from y > to y < with a sign change, 
must have a jump at y = 0. The value of the jump 
$(y = +0) - $(y = -0) = -A/\ui\a xx , and thus the 
boundary values $(y = ±0) = =pA/2\uj\cr xx . 

The formal solution of Eq.(|85|) can be written in 
Fourier components: 



= U\k)(iqA + i X „(k,q)] 
\u>\ [\u>\ + a xx k 2 U(k)] ' 



where k = (k, q), and 
U(k) 



2ire 2 



e(q 2 + k 2 ) 1 / 2 k 



(86) 



(87) 



The constant A is determined from the boundary condi- 
tion: 



dq 
iq J 2tt 



= 0, 



dy^(y 0) = /"ig ( * u (k) - 

7 V \u\<J xx l 

where the second term in the integral is inserted to cancel 
the jump of <!> at y = 0. 

Substituting <f> from Eq. (|8^) , evaluating the part of the 
integral ( |S8| ) containing x^(fc, q) in the limit yo —> 0, and 
simplifying the other part, one obtains 



Xu:(k,q) dq 
a 2tt 



= -A 



\u\ +<J X xU{k)k 2 dq 
\u\ + a xx U(k)(k 2 +q 2 ) 2ir 

(89) 



Now, note that the LHS of Eq.(p9|) is equal to 
* / Xw,k(y)dy = iXuj.k, the one-dimensional source den- 
sity, and the value of $ at y — > is just given by 
— A/2\u>\cr XX) as discussed above. Hence, it follows from 
Eq.® that 



Q(w,fc) 



\uj\ + a xx U{k)k 2 



dq 



<j xx U(k)(k 2 +q 2 ) 2n 



In the special case when Uik) is a constant, the result 
(ppf) agr ees with the expression (59) for q(u>, k) found in 

sec. lira 



The integral over q in Eq.(|90|) for U of the form (pj), 
(|S7| ) can be evaluated exactly. We will only need the 
result for small \k\ -C rj 1 , where r s = e/2irKQ is the 
screening radius of the 2DEG. In this limit, 



TT 



In 



r.|fc| 



+ (1 - a 2 )F(a) 



(91) 



where a = uie/2i:a xx k 1 and 



, (1 — a 2 ) x l 2 arctanVa -2 — 1 



for 
for 



a < 1 
a > 1 

(92) 



The expression (|92j) has no singularity at a = 1. The 
behavior of F(a) as a function of a is such that F(a <C 
1) = tt/2, F(a > 1) = a" 1 La2a, F(l) = 1. 

The next step is to substitute this expression in Eq. ( |79| ) 
to determine the spectral weight AioS) and the instan- 
ton action. The resulting tunneling exponent a(ui) = 
A(uj) + 1 has a weak frequency dependence. This is 
demonstrated on Fig. 0, where a is plotted as a func- 
tion of frequency uj for v = 1/2. In the two other models 
( |S3| ) and (|I]), discussed below, the frequency dependence 
of a{ui) is somewhat stronger. This is quite natural be- 
cause in the model V — V the interaction is to some 
extent screened by image charges, and the results are ex- 
pected to be closer to those for short-range interaction, 
where a(u>) has no frequency dependence. Similar differ- 
ence between the effect of screened and unscreened inter- 
actions on tunneling is known for the diffusive zero-bias 
anomaly MM. 
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FIG. 3. The tunneling exponent a(uj) = A(uA + 1 for the 
models V-V, V + V, and V [see Eqs. (§§), (§j|), and ©] 
at v = 1/2 as a function of frequency to. The frequency is 
measured in units of uo = ftoe 4 - For the model V — V the 
frequency dependence of a is much weaker than for the mod- 
els V + V' and Vo. Note that even in the latter two cases the 
frequency dependence is quite weak, logarithmic at most. 



The model V — V' is closer to the experimental situa- 
tion than other models studied in this paper, because it 
treats interaction as long ranged, and accounts for screen- 
ing in the doped region. Thus, it is the V — V model 
that is interesting to compare to experiment P,p|p9| . The 
tunneling exponent calculated above can be plotted ver- 
sus p xy (see Fig. |^). Experimentally, the parameter con- 
trolling occupation of the Landau levels is the magnetic 
field, and so the experimentally measured a are shown in 
P,|6|Jl9|] as functions of i / l 7 u 1 lk = -B/^o^deg- However, at 
large Hall angle, p xx <C p xy , and away from incompress- 
ible densities, e 2 p xy /h is quite close to v . 

Also, it would be incorrect to ignore the difference 
between the 2DEG density in the bulk and near the 
edge, and to compare the graph in Fig. [| directly with 
the experimentally measured a. One can argue (see 
Sec. [v] below) that the density near the edge exceeds 
I'buik by 20 — 30 %. Taking this into account, one has to 
rescale the slope of the experimentally observed depen- 
dence a = 1/^Wkj and to compare the curves in Fig. || 
with the dependence a — (1.2 — 1.3) p xy e 2 /h. This agrees 
reasonably well with the average slope of the curves in 
Fi g. H in the interval 1 < p xy < 4 studied experimentally 

Of course, a more important issue is whether there 
are plateau-like features in the experimental dependence 
a{v). In the experiment Q a straight line is observed, 
without any sign of plateaus. More recently, however, it 
was found that some samples show signs of a plateau near 
v = 1/3. Upon rescaling of the filling factor by 1.2 — 1.3, 



this corresponds to v e dge between 1/2 and 1/3, which is 
exactly where the middle of the plateau in Fig. ||| is lo- 
cated. However, the matter is clearly not yet resolved, 
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FIG. 4. The tunneling exponent for the model V — V 
[Eq.(js2|)] corresponding to Coulomb interaction screened by 
the doped region in the overgrown edge system. 



There is one other type of interaction for which the 
problem ( |76| ) in the half plane is tractable by Fourier 
transform. It corresponds to the model V + V above, 
defined by Eq. (|83|) . The interaction (83) describes the 
situation when image charges are of the same sign as the 
source charges. Despite being unphysical, this problem 
is still worth attention, because the solution is very sim- 
ple and has the behavior qualitatively different from the 
model V — V'. Physically, this problem is similar to the 
one of unscreened interaction which we discuss below. 



Starting with the interaction (83), one can extend the 
problem to the full plane, now in a symmetric way: 
&(x,—ij) = $(x,y), etc. Upon doing this the interac- 



tion (83) has to become unscreened, of the form U(k) 
given by Eq.(^j). Then the solution is straightforward in 
Fourier components: 



Mk) 



2iXuj,k 



\a xx k 2 +ujyU(k) 



(93) 



This form automatically satisfies the boundary condition 
d y <&(y = 0) = 0, because $ is an even function of y. 

The function Q{uj,k) is found by evaluating $ at the 
boundary y = 0: 



Q- 1 (cj,k) = ± J 



dq 



q 2 + fc 2 + \cj\/a xx U{k) 



irk 



(94) 



where F(a) is defined by Eq.(B3). 
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Again, we now substitute this expression in Eq.(|7E 
to calculate the instanton action. The resulting tunnel- 
ing exponent a(oj) = A(u>) + 1 has a logarithmic fre- 
quency dependence, as shown in Fig. ||. The origin of 
this logarithmic dependence can be traced to the zero- 
bias anomaly in a diffusive conductor p6| , po| . On Fig. || 
we plot a as a function of p xy for several values of lu. One 
notes that the values a in Fig. |^ are somewhat larger than 
those for the model V — V' in Fig. |]. This is due to the 
"antiscreening" in the model V + V which enhances the 
effect of the long-range part of the interaction in the dy- 
namics. Qualitatively, the behavior of a for the model 



C. Wiener-Hopf problem for the model Vo 



Here we consider the model Vq, describing unscreened 
Coulomb interaction (|84|), i.e., in the absence of image 
charges of any kind. The strategy will be to derive an 
integral equation for <&uj,k(y) and to deal with it using 
the Wiener-Hopf method. Our approach is similar to 
that employed by Volkov and Mikhailov in a study of the 
edge magnetoplasmons |^7j . 

We start with the problem (|7^) written in Fourier rep- 
resentation with respect to x. Nondimcnsionalizcd, the 
first equation of (76) reads: 
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FIG. 5. The tunneling exponent for the model V + V' 
[Eq. (^)] corresponding to Coulomb interaction in the pres- 
ence of an "antiscreening" due to image charges in the doped 
region of the same sign as the source charges. 




(fc 2 - <9 2 )<Iv fc (?/) + im u ,k{y) = ix^,kS{y - y ), (95) 

where fi = \w\/cr xx and x = x/(\ UJ \ cr xx)- As m the above 
discussion of the model V — V , it is convenient to place 
the source x(u) aTi a small distance yo from the edge, and 
take the limit yo — > later. 

Posing the correct boundary condition for Eq. (^5|) re- 
quires a discussion. The absence of normal current at 
the edge means that d v Q(y = 0) = 0. On the other 
hand, by integrating Eq.(|95|) from the edge to the source 
X, over the small interval < y < yo, from current con- 
servation one obtains d y <&(y = yo + 0)y a ^o = ~iXu>,k- 
Therefore, in the limit yo — > the boundary condition is 
written as d y $(y — > 0) = — iXui.u- This condition defines 
completely the boundary value problem in the region of 
interest y > yo — > 0. However, without any loss of gener- 
ality, it will be convenient to assume that near the very 
edge, for < y < yo, the normal derivative d y $ vanishes. 

Now, by performing convolution of Eq.(|95|) with 
Uk(y) = J e lkx U(x,y)dx, remembering that nuyy < 0) — 
0, and using the second equation of (76), we transform 
the problem to 



f U k (y-y')(k 2 ~d^. k (y')dy' + ^u J .k(y)=^,kU k (y-yo) (96) 

Jy>>0 

We will be solving Eq. (|9^) in the domain y > with k and u> being parameters. Hence, for simplicity, below we 
suppress the dependence on u> and k and use $(y), U(y), etc. 

It is convenient to integrate in Eq.([9rj) by parts using the boundary condition d y $(y — > 0)j,<j, o = 0, which gives: 

(k 2 -d 2 y ) f U(y-y')$(y')dy' + ti$(y) = ixU(y-y )-d y U(y)$ , (97) 

Jy'>0 



where $o = <&(y = 0). The form (|97| ) of the problem is 
most suitable for applying the Wiener-Hopf method to 
which we now proceed. 

The first step is to perform Fourier expansion of $(y) 
with respect to the y coordinate: 



= (98) 

Since the integral in Eq.(p7|) is taken over y' > 0, in or- 
der to rewrite it in terms of &(q) we decompose $(y) as 
$(y) = $>(y) + $<(y), nonzero for y > and y < 0, re- 
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spectively. One can assume that $>(&/) and $<(y) decay In terms of $> and $<, Eq. ( |96| ) becomes 
at y — > ±00 and verify it later, when solution is found. 

I 

/*(*>(«) + *<(?)) + (k 2 + g 2 )C/(g)$>( g ) = iUfaXxe-*™ - g*o). (99) 



Here the Fourier transformed interaction [/(g) is given 
by (|S7l ). In what follows we set yo = 0. 

The functions ^(g) and <&<(<?) have nice analytical 
properties, namely, "^(g) is an analytic function of q in 
the upper complex half plane Img > 0, and $>(g) is 
analytic in the lower half plane Im q < 0. To make the 
discussion below more transparent, we denote < ! > >(q) by 
$_(g), and $<(g) by $+(g), where ± indicate the half 
plane of analyticity in q. 

Now, Eq. (K)9]) can be written as 



where 



K(q)$_(q) + $ + (q)=R(q), 



K(q) = l + -(k 2 + q 2 )U(q) 
R(q) = -U{q) (X -?*«,). 



(100) 

(101) 
(102) 



The next step is to decompose -ftT(g) into the ratio of two 
functions which are analytic in the upper and lower half 
planes, respectively: 



K{q) = 



where 



X±(q) = exp 



2ni 



X+(q) 

x-{ q y 

dq' 



-00 1 



(103) 



lnK(q')), (104) 



The asymptotic behavior of X±(q) at \q\ ^S> 2ttkq/6 is 
X+(q) = (q + i\k\)/X, X_{q) = X/(q - i\k\), where 
A = ^/JIkq. 

Now, Eq. d9§) turns into 



where 



X+(g) X_(g) 



#(<?) _ «(x-9$o) 



*(?), 



\X.(q) X + {q) 



(105) 



(106) 



Now we decompose ^{q) into the sum of two functions 
with appropriate analytical properties: 



= *+(«) -tf-(g) ; #±(g) = 



1 

2~7ri 



rig' 



The standard Wiener-Hopf reasoning ]28| then leads to 

$ + (< Z )=X+(g)*+(g) ; 



g' - g =F *0 



-X_(g)*_(g). 



*(g). 



(107) 



(108) 



Fourier transform of Eq.( 108 ) gives $(y) for y < and y > 0. 

It is not difficult to find v I / ±(g) explicitly. For that, one has to substitute Eq.( |l06| ) into the Cauchy integral in 
Eq.(107), which gives 



•-<«>— s? 



1 



dq' i( X - g$o) 
-g + iO q 2 + k 2 \X_{q) X+(g) 



(109) 



and a similar equation for \I/ + (g). Now, we close the integration contour in Eq. (]106|) in the upper or lower half plane, 
depending on whether XI 1 or X Z 1 is to be integrated, and evaluate the integral ( pL09[ ) using residues. Having found 
\p_(g), and then using Eq.(|l08|) to go back to $_(g), we obtain 



$ ( ) = *(x~ q *o) X_M 
~ W q 2 + k 2 2\k\ 



1 x + i|*|4 



1 X -i\ k \®o 



i\k\+q X-(-i\k\) i\k\-q X+(i\k\) 



(110) 
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Several remarks are in order about the result ( |11C 

First of all, let us verify tha t Q-(q) is analytic at 
0. The expression (110) has an apparent pole 
-i\k\. However, it is easy 



Imq < 

in the lower hal f pla ne at q = 
to see from Eq.(llO) that the residue for this pole is zero. 
From analyticity at Im q < it follows that $(?/ < 0) = 0, 
as it should be. 

Next, let us verify that the boundary value $o is repro- 
duced correctly. For that we expand Eq.(llO) in inverse 
powers of q at \q\ — > oo: 



*-(«) = - ~ 

iq q A 



(111) 



Since <&(?/ < 0) = 0, one simply has ® (y — > +0) = a. To 
evaluate a, only the first term of Eq.( |110| ) is important, 
because X^(q — > oo) = X/q + Q {q 2 ) ; where A = ^/JTkq, 
and thus the second term of Eq. (|llO| ) does not contribute 




FIG. 6. The tunneling exponent for the model Vo (|84| ) cor- 
responding to unscreened Coulomb interaction. Frequency is 
given in units of ojq = Kge 4 . 



After these consistency checks we can proceed with 
finding the relation between <E>o and x- Conservation 
of current at the boundary y — for the problem 
( |95| ) implies d y <&{y — ► 0) = — i\. On the other hand, 
b = d y Q(y — » 0) in the expa nsion (111). By carrying out 



the expansion of the result (110) up to the order q to 
obtain 6, and then setting up the equation — i\ = b, we 
have 



1 



-IX- 



2\\k\ 



X— i|fc|jo _ x + z|fc|$o 
X + (i\k\) X_(-i\k\) 



, (H2) 



where A is the coefficient in the asymptotic expansion of 
X-(q — > oo) defined above. This equation can be rewrit- 
ten in the form 



_ ix *+(#!)- *-(-#!) 

\k\x + (i\k\)+x_(-i\k\y 



(113) 



According to Eq. (|77|) , the relation (113) defines Q(uj,k) 
in terms of X + (i|/c|) and A_(— i\k\). 

The expressions for X±(±j|A;|) can be simplified: 



X±(±i\k\) = exp [±/(w, fc)] 



(114) 



where 



I(uJ,k) 



in (i + -k\e + i)u{k^fe+i\ 

(115) 



Here £ = ^ = /cr^ 

(113), one finally arrives at 



After putting (|114|) into 



Q(w, fc) = |fc|cothJ(w,fc) 



(116) 



With this expression for Q(ui, k) one can go back to the 
effective action (|78|), and find the Green's function (|79| ) 
in terms of the spectral wei ght A (ui) given by (|80|). 

The integral entering Eq. (|ll5| ) can easily be tabulated 
numerically. The spectral weight A(uj) has a logarithmic 
frequency dependence, as shown in Fig. ||[ similar to that 
of the model V + V . The behavior of the tunneling ex- 
ponent a as a function of p xy , shown in Fig. |(], is also 
close to that for the model V + V. One notes that the 
values of a are somewhat less than those for the model 
V+V with similar parameters. This is due to a relatively 
weaker effect of the long-range part of the interaction in 
the model Vq. 



V. COMPARISON TO THE EXPERIMENT 

In this section we discuss some aspects of the over- 
grown cleaved edge system J|,||,|l9|]. In our view, the 
most relevant issue concerns the 2DEG density distri- 
bution near the edge. One of the key features of cleaved 
edge systems is that they produce structures with sup- 
posedly an atomically sharp confining potential, and thus 
the 2DEG density profile near the edge is expected to be 
reasonably smooth. This is important in edge tunneling 
experiment, because the system must have a well defined 
filling factor even very close to the edge. 



A. Thomas-Fermi model 

To estimate the importance of various factors control- 
ling the density near the edge, below we consider a sim- 
plified electrostatic Thomas-Fermi model, in which the 
2DEG is modeled as an ideal charge fluid, and all effects 
of electron-electron correlation and finite density of states 
are ignored, except very close to the edge. In principle, 
this approximation is quite reliable at distances larger 
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than the screening length r s = e/2-kkq, and so the re- 
sults will be meaningful at distances more than r s from 
the edge. 

The electrostatic problem we consider involves the 
2DEG density n(x, y) in the half plane y > 0, top surface 
charge states that are at a distance w = 600 nm above the 
2DEG, a layer of charged donors parallel to the 2DEG 
at a distance w + = 60 nm above the 2DEG plane, and 
also charges in the three-dimensional doped region, which 
in our model occupies the halfspace y < —Wb, where 
W(, = 9 nm is the width of the barrier together with the 
buffer region. The top surface, the 2DEG, and the doped 
region are assumed to be equipotentials in the problem. 
For simplicity, we assume that the 2DEG is grounded, 
and the bias voltage on the 3D doped region is very small, 
so that the electrochemical potentials of the two regions 
are essentially equal. Relative to the 2DEG, the elec- 
trostatic potential at the top surface is V s — — 800 mV, 
and the electrostatic potential at the boundary of the 
3D doped region is Vd ~ 20 mV. (The value of Vd re- 
flects the chemical potential difference before the charge 
redistributes itself. It is given by the difference of Fermi 
energies in the doped region and in the 2DEG plus the 
confinement energy of the 2DEG.) The charge density of 
donors er+ is taken to be constant everywhere at y > up 
to the edge y = 0. The potential Vd is much smaller than 
the barrier height, which is estimated as ~ 120 meV. 

One can write down a simple analytic formula for the 
2DEG density, using the electrostatic superposition prin- 
ciple, according to which the effects on the 2DEG due to 
the donors, the top surface charge, and the doped region, 
can be treated separately and then added. 

First, let us consider the charge induced by donors, 
when the top surface and the doped region are at the 
same electrostatic potential as the 2DEG. We make an 
approximation w 3> u>+, which allows us to move the top 
surface to infinity, and thus to ignore it. Also, we assume 
that the distance to the doped region Wb <C w + , the sep- 
aration of the donors from the 2DEG. With the values 
for w, W-\-, and Wb quoted above, both approximations 
are reasonable. The resulting contribution to the 2DEG 
charge density is 



7 2DE G (y) = ^ arctan 



(117) 



It describes the 2DEG density, constant and equal to a + 
at y ^> w + , and decreasing to near the edge. 

The effect of the top surface potential V s , in the ab- 
sence of donors, and with the 2DEG and the doped re- 
gion at zero electrostatic potential, can be evaluated as 
follows. In the approximation Wb <C w, the problem is 
equivalent to the standard electrostatic problem of a half- 
open slit, with one side of the slit being at the potential 
V s with respect to the other side and the end. The in- 
duced charge density in this problem is 



2DEG 



(y) 



P- tanh H. 

Anw 2w 



(118) 



This contribution is constant and equal to cV s /4ttw in 
the bulk, at y 3> w/tt, and decreases to zero near the 
edge. 

Finally, the effect of potential difference between the 
2DEG and the doped region can be considered ignoring 
the top surface and the donors. The relevant spatial scale 
in this case is Wb <C w + , w, and so the problem is reduced 
to that of a ground half plane (representing the 2DEG), 
and a conducting plane perpendicular to it, at a relative 
potential Vd, located a distance Wb away from the ground 
half plane. The charge density induced in the 2DEG is 



1 



Wb) 



(119) 



It behaves as 1/y away from the edge, and as lf-y/y near 
the edge. The square root divergence near the edge is 
an artifact of the simplified model ignoring finite den- 
sity of states of the 2DEG. In a Thomas-Fermi model, 
the divergence would be cut at a distance ~ r s from the 
edge. 




100 150 200 250 300 350 400 450 
distance from the edge, nm 

FIG. 7. Density distribution in the 2DEG near the edge 
plotted for six values of the potential Vd of the doped re- 
gion (listed from top to bottom). The top surface potential 
V s = -800 mV and the donor density a+ = 1.94 • 10" 
correspond to the 2DEG bulk density ribuik = 10 11 cm' 
geometrical parameters used are defined in the text. 
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The resulting 2DEG charge density is a sum of three 
terms, cr tatal = cr^ EG + ct^deg + °2DEG ■ To eliminate 
the unphysical singularity near the edge due to cr' 3 ' , we 
average the density Ototai over intervals of length 2r s , and 
consider 



'total 



(V) = (2r.)" 



(Ttotai(y')dy' . 



(120) 
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The averaged density crj 1 ™^ is plotted in Fig. [?] for several 
values of the doped region potential Vd- The screening 
radius used in the averaging is taken to be r s — 20 nm. 

One can see from Fig. [7] that the density within ~ 
200 nm near the edge is quite sensitive to the potential 
Vd- Another feature evident in Fig. [?] is that the density 
close to the edge exceeds that in the bulk by 20 — 30 %. 
The 2DEG density approaches the bulk value at distances 
> 400 nm from the edge. Also, there is a peak in the den- 
sity profile near the very edge, resulting from the 
contribution averaged over the length ~ r s . This peak 
makes the density profile nonmonotonic, with a minimum 
at ~ 30 - 40 nm from the edge. Altogether, the 2DEG 
density near the edge is smooth but not perfectly uni- 
form. 

It should be remarked that our simplified electrostatic 
model is perhaps insufficient at distances smaller than or 
of the order of r s ~ 20 nm. Thus the smallest scale fea- 
tures in Fig. 0, such as the density peak near the edge, 
should be taken with caution. Moreover, we used the 
Thomas- Fermi model, the screening radius r s , etc., in the 
absence of magnetic field. It remains to be seen whether 
the results are preserved in a more accurate treatment 
accounting for Landau levels, finite 2DEG compressibil- 
ity, and exchange effects. On the other hand, on spatial 
scales larger than r s , the results obtained from a purely 
electrostatic model should be reliable. 

One issue that can be addressed using the electrostatic 
model is the calibration of density in the experiment |^|,^| . 
The tunneling exponent a is presented there as a func- 
tion of magnetic field, which is calibrated in terms of the 
bulk filling factor using magnetotransport data. How- 
ever, the filling factor relevant for tunneling is that near 
the edge. According to the above, in the region 100 ~ 300 
nm from the edge, the density is at least 20 — 30% higher 
than in the bulk. If one assumes that this is the relevant 
distance scale for charge relaxation at the temperatures 
and voltages employed in the experiments, then the de- 
pendence a — 1/j/buUc observed in [|6|Jl9|l translates into 
a w (1.2 — 1.3) z^jgc- In actuality, the relevant distance 
scale will depend on the filling factor and the cleanliness 

I 




of the edge, as well as the energy of the tunneling elec- 
tron. 

One notes that after accounting for the difference be- 
tween i/ e dge and t'buik the dependence a{v) shifts closer 
to the theoretical curves (see Fig. ||). 



B. Two-mode model 

Because the 2DEG density profile discussed above is 
significantly nonmonotonic near the edge, it is possible 
that this may change the structure of the edge modes. 
More precisely, suppose that the peak density near the 
edge is so high that the filling factor reaches v = 1 within 
the region ~ 30 nm corresponding to the peak displayed 
in Fig. f^. Then the edge mode on the periphery will 
correspond to v = 1 even when v < 1 away from the 
edge. In this case, in addition, there will also be coun- 
terpropagating modes positioned on the inner side of the 
incompressible v = 1 region. The number of these modes 
and their Hamiltonian will depend on v somewhat away 
from the edge. This type of acomposite structure of the 
edge was first proposed by MacDonald for the v = 2/3 
system, based on a Hartree-Fock analysis |^9| . 

In this model, the tunneling electron is injected into 
the outer v = 1 mode, because of higher overlap of the 
tunneling state with the mode closest to the edge. We 
assume that the edge is so clean that we can neglect scat- 
tering between different edge modes. Then, the inner 
modes will be important only to the extent that tunnel- 
ing charge couples with them by Coulomb interaction, 
and shakes them up. In this scenario, after tunneling 
there is no statistics change of the injected particle, since 
it remains in the fermionic v = 1 edge state. There- 
fore, one expects a smooth dependence of the tunneling 
exponent on v, without any cusps or plateaus. 

To estimate the shakeup effect due to Coulomb cou- 
pling to the inner modes, let us represent them by a sin- 
gle charged mode. Thus the system can be described by 
two counterpropagating chiral modes: 



iujk 



?< /> (2 V ( 



2) 



(121) 



where g = 1 — v and Vy is the coupling matrix, expressed 
in terms of the interactions VR as follows: V\\ = V® x , 
Via = g ffi, V 2 i = gV§ x , V 22 = g 2 V 2 ° 2 . The form of the 
action (|l2l|) can easily be justified in the same way as 
in Sec. |III|. In this case there is no issue of charge injec- 
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G(r) 



exp 



1 

8^ 



tion in the inner mode, and so there are no complications 
related with counterterms, as in (|l2"|). 

It is straightforward to write down the Green's func- 
tion by evaluating the saddle point of the quadratic ac- 
tion (121). The result reads 



(tj 2 V22 — igcjk)\Ju\ 2 did dk 
(w 2 Vii + iojk)(uj 2 V22 — igujk) — cj 4 Vi2V2i 



(122) 



23 



To evaluate the Green's function, we assume that the 
coupling matrix Vij has no k dependence. This is true for 
the screened Coulomb interaction 2-7re 2 (l — e~ 2a ' fe ')/e|&;| 
at a\k\ < 1, where a is the distance from the edge mode 
location to the doped region. Hence the length a is some- 
what larger than the barrier width Wb- 

In this case the integral over k can be done by residues, 
and the result is G(t) = r~ Q , where 



V? 
v ll 



gv 2 ° 2 



(123) 



The dependence a(u) in the interval < v < 1 is smooth, 
without singularities, as it should be in the case when the 
effect of the fractional edge is purely a shakeup, not ac- 
companied by injection of charge. 

To estimate numerical values of a, we consider a model 
in which the interactions are given by the Coulomb 
potential screened by the doped region. We assume that 
the outer edge is separated from the doped region by 
a barrier of thickness u>, and the inner and outer edge 
states are a distance a apart. Then 
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We consider the limit of small fc, where the interac- 



tions V^ ' do not depend on k: V^"' = 2w(2ite 2 /e) 



(o) _ 



V$ = 2(w+a)(27re 2 /e), V 1 



(0) 
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(2w-a)(2ne 2 /e). 
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FIG. 



Solid curv es: The tunneling exponent in the 



composite edge model (121) is shown for three values of the 



ratio a/w of the distance between the outer and inner edges 
and the tunneling barrier width. For comparison, a theoret- 
ical curve for the V — V' model is shown, for p xx = p xy /W 
and uj = 10 _5 Koe 4 . Dashed curves: The straight line 
ct = p xy corresponds to experiment (Refs. [^|],|l!|), and the 
line a = 1.3p xy is obtained by correcting the filling factor by 
the ratio of the densities near the edge and in the bulk. 

In this model, the only parameter is the ratio a/w. 
The tunneling exponent is plotted in Fig. as a function 



of p xy = 



for several values of a/w. On the same 



I 

figure, we show the experimental dependence of a versus 
p xy rescaled by a factor 1.3 as discussed above. 

The distinct feature of the composite edge model is 
the absence of plateaus in the tunneling exponent a{v). 
However, note that in order for the tunneling exponent 
a to fall in the right range, one has to assume unphysi- 
cally small values of the ratio a/w. Also, the theoretical 
curves for nonzero a/w have curvature which is absent 
in the experimental curve. This curvature is even more 
significant at higher values of the parameter a/w and is 
unlikely to disappear if one takes into account possible 
dependence of a/w on v. It is apparent that this sim- 
plified two-mode model does not agree with the experi- 
mental results on tunneling. Nevertheless, it illustrates 
the point that, if scattering between edge modes is suf- 
ficiently small, a complicated edge structure can lead to 
large changes in the observed tunneling exponent, which 
will not be closely related to the bulk filling factor. 



VI. SUMMARY 

The problem of tunneling into the edge of a composite 
fcrmion QH system is treated for long-range Coulomb in- 
teraction between electrons, as well as for a short-range 
interaction model. It is shown that in the case of diffu- 
sive CF dynamics described by finite p xx , the tunneling 
exponent is controlled by the coupling of tunneling elec- 
tron to the charged edge mode. The effective action for 
this mode is a generalized chiral Luttinger action with a 
nonlocal dissipative term. 

The tunneling exponent is found to be a continu- 
ous and monotonic function of p xy , given, in the limit 

Px X — » 0, by a — 1 + ^(p xy - \pxy\), where p x °y is the CF 
Hall resistivity due to motion in the residual magnetic 
field. In order to verify the robustness of the results we 



24 



consider several models for the electron-electron interac- 
tion: the short-range and Coulomb interaction, and, in 
the latter case, with and without electrostatic screening 
due to image charges in the doped region. 

The dependence of a on p xy is characterized by 
plateaulike features, not observed in the experiments on 
cleaved edge systems. We discuss the 2DEG density pro- 
file near the cleaved edge, and propose that the discrep- 
ancy between theory and experiment is possibly due to 
spatial variation of the density near the edge and, in par- 
ticular, to a nonmonotonic density profile, giving rise to 
a composite structure of edge modes. 
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APPENDIX A: ESTIMATE FOR THE BALLISTIC 
REGIME 

In the ballistic regime, for length scales smaller than 
the composite fermion mean free path I, the conductiv- 
ity tensor is nonlocal in space. Close to the edge, the CF 
conductivity er^g (r, r', ui) depends on the distance from 
the edge, as well as on the separation r — r'. As a crude 
approximation, however, in order to estimate the con- 
tribution of the short-distance response to the tunneling 
exponent, we shall ignore the dependence on the distance 
from the edge and use, instead, the bulk CF conductivity, 
which depends only on r - r' |25[|. 

As discussed in Sees. || and ID, for a nonlocal con- 
ductivity we may still approximate the Green's function 
G(r) using the factorization (|63|), but the actions S(t) 
and S l f rcc (r) should be evaluated using the correct nonlo- 
cal conductivity. Instead of this, in our approximation, 
we use the form ( |65| ) for S, with the change that we re- 
place the macroscopic conductivity a xx by the quantity 
<7;(|fc|), which is the wave- vector dependent longitudinal 
conductivity for the bulk compressible Hall state. Specif- 
ically, at v = 1/2, according to Ref. |Q, we have 



tnik) ~ (e 2 /8Trh)(k/k F ), 



(Al) 



for l^ 1 < k < k jp , while er;(fc) reduces to the macroscopic 
conductivity a xx for k < l^ 1 . We continue to approxi- 
mate the Hall conductivity in Eq. ( |65| ) by its macroscopic 
value a xy . 

In evaluating Eq. (|6^) , it is convenient to combine con- 
tributions from wave vectors k and —k, and replace the 
sum over k by an integral over positive values of k. If 



the frequency u> is sufficiently small, there will be two 
distinct regions which can contribute significantly to the 
integral. The first, for k < £ , gives the same contribu- 
tion as was found in Sec. Ill above, since <Ji(k) = a xx 



in this region. In the region k > £ , we may set q = fc, 
when lo is sufficiently small, so that the integrand takes 
the form 



h 



largck 



\J{w)f 



-Re 



Or 



\u\\k\(j xy (<Ji(k) + ia xy ) ' 



(A2) 



The last factor in Eq.(A2) is small for k <^_kp but be- 
comes of order unity for k ss k F , where ui{k) « a xy . The 
contribution from this region to the integral could there- 
fore make a contribution of order unity to the tunneling 
exponent a. However, this contribution may be largely 
or completely canceled by the corresponding contribution 
to S[ TCC . 

If we neglect the difference between the longitudinal 
and transverse conductivities at the finite wave vector fc, 
then 



Re 



1 



(cr/(fc) + ia 



Pi{k) 



(A3) 



xy) 



which is the longitudinal resistivity at wave vector k. In 
calculating S[ rcc , using the same assumptions, we obtain 
the identical expression, because the longitudinal conduc- 
tivity of the composite fermions is the same as that of the 
electrons. Thus the contribution to the tunneling expo- 
nent from short wavelengths is canceled, in this approx- 
imation. We therefore wind up with the same value for 
a as was obtained in Sec. Ill, namely a ~ 3 at v = 1/2, 
for a system where a xx <C a xy . 

It is not possible to say whether a similar cancelation 
would occur in a proper analysis incorporating the nonlo- 
cal conductivity. If the cancellation does not occur, then 
the surviving contribution from short wavelengths could 
give a contribution of order unity to the tunneling expo- 
nent, which would be independent of the mean- free-path 
in the limit I <C kp. 
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